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ANALYTICAL STUDY OF EFFECTS OF PRODUCT OF INERTIA 


ON AIRPLANE SPIN ENTRIES, DEVELOPED SPINS, 

AND SPIN RECOVERIES 

By Ernie L. Anglin 
Langley Research Center 


SUMMARY 


An analytical study has heen made to investigate the effects of variations 
in product of inertia on the spin entry, developed spin, and spin recovery of 
airplanes. Three configurations considered to he representative of modern 
fighter-type airplanes were used: a delta-wing fighter, a swept-wing fighter, 

and a stub-wing research vehicle. They were mass loaded relatively heavily 
along the fuselage. 

The results indicate that the magnitude of the product of inertia may have 
a significant effect on attempted spin-entry motions. For the delta-wing 
fighter and the swept-wing fighter, an increase in product of inertia changed 
the attempted entry motion from one in which a developed spin was attained to 
one in which no spin resulted. For developed spins, increasing the product of 
inertia caused slightly higher angles of attack, somewhat lower rates of rota- 
tion, and lower magnitudes of spin-energy factor resulting in fewer turns 
required for recovery. When spins are oscillatory, the trends noted for devel- 
oped spins are somewhat obscured. 


INTRODUCTION 


It has long been recognized (for example, see refs. 1 and 2) that moment- 
of-inertia variations have large effects on developed- spin and recovery motions 
of airplanes and that these motions are particularly affected by whether the 
vehicle is mass loaded more heavily along its fuselage or its wing and by the 
degree of difference between the wing and fuselage loading. Results of refer- 
ence 1 show that this difference in wing and fuselage loadings is extremely 
important both in algebraic sign and in magnitude. No information exists, how- 
ever, with regard to possible effects of product-of -inertia variations due to 
varying displacements of the principal longitudinal axis with respect to the 
fuselage reference (body) axis of the airplane. 

The present investigation was therefore made to determine analytically the 
effects of product -of -inertia variations on the spin-entry, developed- spin, and 
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spin-recovery motions. The product -of -inertia variations used correspond to 
angular displacements between the reference and principal axes from 0° to 5°« 
Three representative modem airplane configurations mass loaded relatively 
heavily along the fuselage were used: a delta-wing fighter, a swept-wing 

fighter, and a stub-wing research vehicle. A high-speed digital computer, 
wind-tunnel aerodynamic data, six-degree-of-freedom equations of motion, and 
methods similar to those used in references 3 , 4, and 5 were utilized to cal- 
culate the airplane motions. 


SYMBOLS 


The body system of axes was used in the calculations. This system of axes, 
related angles, and positive directions of corresponding forces and moments are 
illustrated in figure 1. 
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side-force coefficient due to aileron deflection, per deg 
side-force coefficient due to rudder deflection, per deg 
normal-force coefficient due to elevator deflection, per deg 
mean aerodynamic chord, ft 

i • 2 

\ 

non dimensional spin energy factor, - 

I pV R 2 Sb 

longitudinal force acting along X-axis, lb 
side force acting along Y-axis, lb 
normal force acting along Z-axis, lb 

acceleration due to gravity, 32.17 ft/sec^ 
altitude, ft 
initial altitude , ft 

moment of inertia about vertical axis, slug-ft^ 

moment of inertia about Xq-, Yq-, and Z Q -axis, respectively, 
slug-ft 2 

moment of inertia about X-, Y-, and Z-axis, respectively, 
slug-ft 2 

product of inertia in XZ-plane, positive when X 0 -axis is 
inclined below X-axis at nose, slug-ft^ 

inertia yawing-moment parameter 

rolling moment acting about X-axis, ft -lb 
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pitching moment acting about Y-axis, ft-lb 


yawing moment acting about Z-axis, ft-lb 
mass, W/g, slugs 

component of resultant angular velocity about X-, Y-, and 
Z-axis, respectively, rad/sec 

wing area, sq ft 

time, sec 

component of resultant linear velocity Vp along X~, Y-, and 
Z-axis, respectively, ft/sec 

resultant linear velocity, ft/sec 

weight, lb 
body axes 

principal axes (axes about which the products of inertia are 
zero) 

angle of attack, angle between relative wind Vr projected 

into XZ-plane of symmetry and X-axis, positive when relative 
wind comes from below XY -plane, deg 

angle of sideslip, angle between relative wind V R and projec- 
tion of relative wind on XZ-plane, positive when relative 
wind comes from right of plane of symmetry, deg 

aileron deflection with respect to chord line of wing, positive 
when trailing edge of right aileron down (left stick), deg 

elevator deflection with respect to fuselage reference line, 
positive with trailing edge down, deg 

rudder deflection with respect to fin, positive when trailing 
edge left, deg 

angle between X-axis and XQ-axis, positive when principal axis 
is below reference axis at nose, deg 

total angular movement of X-axis from horizontal plane measured 
in vertical plane, positive when airplane nose is above hori- 
zontal plane, rad or deg 
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p air density, slugs/cu ft 

0 angle between Y-axis and horizontal measured in vertical plane, 

positive for erect spins when right wing downward and for 
inverted spins when left wing downward, rad or deg 

0 e total angular movement of Y-axis from horizontal plane measured 

in YZ-plane, positive when clockwise as viewed from rear of 
airplane (if X-axis is vertical, 0 e is measured from a ref- 
erence position in horizontal plane), rad or deg 

horizontal component of angular deflection of X-axis from ref- 
erence position in horizontal plane, positive when clockwise 
as viewed from vertically above airplane, rad or deg 

A dot over a symbol represents a derivative with respect to time; for 
example, u = du/dt. 


METHODS ADD CALCULATIONS 


The spin-entry attempts and the developed- spin and spin-recovery motions 
were calculated by using a high-speed digital computer which solved the equa- 
tions of motion and associated formulas listed in the appendix. The equations 
of motion represent six degrees of freedom along and about the body system of 
axes. (See fig. 1 for illustration of body axes.) 

Sketches of the planforms of the three configurations used in the study 
are shown in figure 2. Configuration A represents a delta-wing fighter, config- 
uration B represents a swept-wing fighter, and configuration C represents a 
stub-wing research vehicle. The mass and dimensional characteristics of these 
configurations are listed in table I. 


The aerodynamic data inputs used in the digital computer calculations are 
presented in figures 3 to 8. The data for configurations A, B, and C were 
obtained from references 4, 6, and 3, respectively. Estimated values of the 
derivatives C^ and C^ were used in the pitching equation of motion (see 

appendix) and were constant for all angles of attack. Values of C-^ and 


Cm<x vere -0.45 and -0.45 for configuration A (ref. 4) and -10.0 and 0 for con- 
figuration C (ref. 3). Values of C^ and C^ used for configuration B in 


this investigation were estimated to be -2.0 and 0, respectively. No values of 
the derivatives Cj and C n for configurations B and C were available and 

for purposes of this study they were assumed to be 0. 


The procedure for calculating attempted spin entries was as follows: The 

configuration initially was in trimmed level flight. At zero time, the eleva- 
tor was deflected full up from its trimmed level flight position. As the motion 
developed, rudder and aileron controls were applied at times conducive to the 
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attaining and sustaining of a developed right spin (turning toward the pilot 1 s 
right). The initial conditions used for these calculations are shown in 
table II. 

Developed spins were obtained in a manner similar to the spin-tunnel 
testing technique. (See ref. k.) In this method, it is assumed that the 
initial condition is one of very high angle of attack with applied rotation 
about a vertical axis and that from this condition the spin parameters will 
undergo some intermediate motions until the configuration eventually achieves 
its own equilibrium-developed spin. The initial conditions used for these cal- 
culations are included in table II. 

Spin-recovery attempts were made by deflecting the rudder against the 
direction of rotation and by deflecting the ailerons with the direction of roll 
(left rudder and right ailerons when in an erect spin to the pilot f s right 
these are the optimum control deflections for recovery from developed spins for 
airplanes loaded relatively heavily along the fuselage (ref. l), as are the con- 
figurations investigated herein. The elevators remained in the original 
elevator-up setting. A spin is considered terminated when either the spin 
rotation ceases or the angle of attack becomes and remains less than the stall 
angle. Usually when the angle of attack becomes less than the stall angle, the 
airplane enters a steep dive without significant rotation (\jr e « 0) . In some 
instances, however, the airplane may be turning or rolling in a spiral glide or 
an aileron roll. Also, sometimes the airplane may roll or pitch to an inverted 
attitude from the erect spin and may still have some rotation, but it is con- 
sidered to be out of the original erect spin. 

The product-of-inertia variations were assumed to be due to weight shifts 
inside the aircraft, accomplished in such a manner as to cause angular displace- 
ments of the XQ-axis (principal axis) below the reference X-axis (body axis) 
at angles ranging from 0° to 5°. The equations used to obtain values for the 
moments and product of inertia about the body axes are from reference 7 and are 
as follows: 

% - §(% + IXo) - !(lz 0 - Ix 0 )c°s 2e 
= IY 0 

iz - |(iz 0 + IXo) + - I X 0 ) COS * 

X XZ = - x x 0 ) sln 2e 


The moments and product of inertia obtained from these equations for each con- 
figuration are presented in table III as a function of e. In addition, values 


of the inertia yawing-moment parameter 


T X ~ *Y 
mb 2 


for each configuration are 


i 
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given in table III. Changes in the algebraic sign or large changes in the mag- 
nitude of this parameter may have an important influence not only on the devel- 
oped spin characteristics but also on the particular combination of controls 
which are considered to be the optimum controls for recovery (ref. 1). The 
variations in yawing-moment parameter shown in table III are not of a large 
enough magnitude to be important. 


RESULTS AND DISCUSSION 


Results of the calculations of the spin-entry, developed-spin, and spin- 
recovery motions for e = 0° and e = 5° are presented in figures 9 to l4 as 
time histories of angle of attack, angle of pitch, sideslip angle, roll atti- 
tude, rate of spin rotation, number of spinning turns made, and aileron- and 
rudder-control inputs. Some pertinent results from the developed-spin and spin- 
recovery calculations are also presented in table IV for values of e from 0° 
to 5° in 1° increments . 


Spin-Entry Calculations 

Spin-entry attempts were calculated for each configuration for values of e 
of 0° and 5 ° only. Therefore, only overall effects were obtained — in partic- 
ular, whether the variation in magnitude of product of inertia used would change 
the attempted entry motion from one in which a developed spin was attained to 
one in which no spin resulted. The same initial conditions and control move- 
ments were used for both the e = 0° and the e = 5° calculations for each 
model . 

Configuration A .- The results of the spin-entry calculations for configura- 
tion A are presented in figure 9 . For e = 0° (ixz = 0 slug-ft 2 ), a spin-entry 
motion occurred in which two turns were made in the first 27 seconds. The angle 
of attack increased smoothly upward with oscillations steadily decreasing in 
amplitude, and the configuration apparently was headed toward an equilibrium 
developed spin condition with an angle of attack above 6o° and with a rate of 
spin rotation \jr e of less than 1 radian per second. 

For e = 5 ° ^XZ = slug-ft 2 ), the angle -of -attack curve shows two 

oscillations which were approximately the same as those obtained for e = 0°. 
After that, however, the angle of attack for e = 5 ° increased very rapidly 
until it went above 120° and beyond the aerodynamic data range used in this 
study. At the same time, 0 e and p were changing violently and ijr e did not 

appear to be headed toward any positive equilibrium value. Although the motion 
did not show beyond question that this is a no-spin condition, past experience 
indicates that an airplane normally would not enter a spin after having gone 
through such gyrations. 

Configuration B . - The results of the spin-entry calculations for configura- 
tion B are presented in figure 10. For e - 0° the motion was oscillatory and 
a rollover occurred after approximately 10 to 11 seconds, after which an erect 
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developed spin was obtained. The angle~of -attack curves for e = 0° and 
e = 5° coincide for the first two oscillations, after which oscillations for 
g = 5° occurred about a much lower angle of attack. Finally, the angle of 
attack became negative; this indicates a no-spin condition. 

Configuration C . - The results of the spin-entry calculations for configura- 
tion C are presented in figure' 11. Spins were obtained with angles of attack 
greater than 8o° for both e - 0° and e = 5°.? and increasing e had little 
or no effect on the entry motions. The spin entry did progress slightly faster 
for € = 0°, however, inasmuch as two turns were achieved at 31 seconds com- 
pared with 34 seconds for e = 5°* 

Summary of entry results .- A variation of 5° th e displacement of the 

XQ-axis (principal axis) below the reference X-axis (body axis) may have a 
significant effect on attempted spin entries. In general, an increase in e 
appeared to increase the degree of the oscillations occurring during the post- 
stall gyrations. Sometimes these increases in oscillations may lessen the tend- 
ency of a given configuration to enter a spin. 

The time histories of configuration C show very mild oscillations when com- 
pared with those obtained for configurations A and B. (Compare fig. 11 with 
figs. 9 a nd 10.) An examination of the equations of motion (see appendix) indi- 
cates that the magnitude of the effect produced by the product of inertia I 
depends on the magnitudes of p, q, and r. A larger effect would therefore 
be expected on motions which are more oscillatory. 


Developed-Spin Calculations 

The results of the developed- spin calculations are presented in table IV 
for configurations A, B, and C for values of e of 0°, 1°, 2°, 3 °, 4°, and 5°* 
The same control settings and initial conditions were used for all values of e 
for each configuration. The calculations for each configuration allowed time 
for any disturbance caused by the specific combinations of initial conditions 
used to damp out and for the configuration to achieve its own equilibrium devel- 
oped spin. The spin parameters presented in table IV are for the equilibrium 
developed spins and are, therefore, the values present just prior to the appli- 
cation of the recovery controls. Recovery controls were then applied. 

Configuration A .- The results of the developed-spin calculations for 
e = 0° and e = 5° for configuration A are presented in figure 12. With 
€ = 0° (fig. 12(a)), configuration A had a very steady developed spin with an 
angle of attack of 74°, sideslip and roll angles near 0°, and a rotation rate of 
1.26 radians per second. The results presented in table IV indicate the manner 
in which the developed-spin parameters are affected as e was varied from 0° 
to 5°. A comparison of figures 12(a) and 12(b) (or the developed- spin results, 
for e = 0° and e - 5° shown in table IV) indicates that the steady nature of 
the developed spins was not affected by changes in e. The angles of attack 
remained at approximately 74°, although table IV shows that a very slight 
increase occurred as € increased, the sideslip and roll angles remained near 
0°; and the rate of rotation decreased from 1.26 for e = 0° to 1.15 for 
€ = 5 °. 
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Configuration B .- The results of the developed-spin calculations for con- 
figuration B are presented in figure 13- With e = 0°, configuration B had an 
oscillatory developed spin. The angle of attack oscillated from 65° to 85° , 
the sideslip and roll angles oscillated from approximately -20° to 20°, and the 
rate of rotation varied from 1 .69 to 2.63 radians per second. 

An examination of the results presented in table IV indicates that as € 
increased, the angle-of-attack oscillations increased, whereas the sideslip- 
angle and roll-angle oscillations decreased a few degrees. The average angle 
of attack remained at approximately 76 °, although table IV shows that a very 
slight increase occurred as e increased. However, the trend was not as con- 
sistent for configuration B as it was for configuration A. The average rate of 
rotation decreased from approximately 2.17 for e = 0° to 2.00 for € = 5°j a 
trend similar to that noted for configuration A. 

Configuration C.- The results of the developed-spin calculations for con- 
figuration C are presented in figure l4. For e = 0°, configuration C had a 
relatively steady developed spin with an angle of attack of approximately 85° 
and a rate of rotation of 2.07 radians per second. 

An examination of table IV indicates that the steady nature of the devel- 
oped spin remained the same for all values of e. The average angle of attack 
remained at approximately 85°, although table IV shows that a very slight 
increase occurred as € increased. The rate of rotation decreased from 2.07 
for € = 0° to I.65 for e = 5°- 

Summary of developed-spin results .- A summary of the results of the 
developed- spin calculations is presented in figures 15 and l6 where the values 
of average anc ^ ^e, average; respectively, are plotted as a function of e. 

The average angle of attack (fig. 15) Tor each configuration remains approxi- 
mately the same as c increases. Very slight angle-of-attack increases did 
occur, but the magnitude of these increases is too small to be significant. 

The rate of rotation (fig. 1 6) for each configuration decreases as € increases. 

Analysis of developed-spin results .- The magnitude of the effect produced 
by product of inertia depends on the magnitudes of p, q, and r. In general, 
in a developed spin, the magnitudes of the oscillations and, therefore, the 
magnitudes of p, q, and r are not as great as those which occur during a 
spin-entry maneuver. Hence, the magnitude of the effect produced by a product 
of inertia would be expected to be less on a developed-spin motion than on a 
spin-entry motion. 

The probable effect of product of inertia on the relative difficulty of 
spin recovery may be determined from the following analysis. Since the average 
angle of attack remained approximately the same when e was increased, the 
decrease in as e was increased may be explained by an equation from ref- 

erence 1. This equation, which approximates the rate of rotation of a vehicle 
in a developed spin, is 
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When this relationship is revised to include the product of inertia, it becomes 



-My 


|( X Z - x x) sln 2a " x xz 


cos 2a 


Inasmuch as the average angle of attack remains approximately the same, the 
numerator of this expression remains approximately the same. However, as € 
increases, the magnitude of the denominator increases for all developed spins 
having an angle of attack greater than 45° (a > 4-5° for each of the configura- 
tions investigated herein); this increase causes a decrease in the value of ty e . 

Since the average angle of attack remains approximately the same as € 
increases, the accompanying decrease in rotation rate becomes important. As 
was pointed out in references 1 and 5* the lower the rate of spin rotation, the 
more rapid the recovery should be. In reference 5* a nondimensional spin-energy 
factor E s , based on the kinetic energy of spin rotation, was shown to indicate 
the relative difficulty of spin recovery. The equation used for calculating E s 
in reference 5 revised for this investigation to include the product of inertia 
becomes 

1 • 2 

2 V^e 

E s = £ 

S 1 P 

I pVr Sb 

where 

Iy - Iy cos 2 a + I £ sin 2 a - sin 2a 


Values of E s were obtained for configurations A, B, and C by using the 
revised equations and are presented as a function of e in figure 17- For all 
three configurations investigated, the magnitude of E s generally decreased 
consistently as € increased. The decreases in E s for configuration A are 
small in terms of absolute values but are large in terms of percent change, and 
actually represent a significant reduction. This reduction is illustrated more 

E 

effectively in figure l8 where the variation of — — - — with € is 

J e E s for € = 0° 

presented. The values of E s for e = 5° represent approximately 8o, 8l, and 

64 percent of the E s values obtained for € = 0° for configurations A, B, 

and C, respectively. 
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Spin-Recovery Calculations 

Attempts were made to calculate the spin-recovery from each of the 
developed- spin calculations. The results of the spin-recovery calculations are 
included in figures 12 to Ik and table IV. 

A summary of the calculated recovery results is presented in figure 19, 
in which turns for recovery are presented as a function of e . The turns 
required for recovery for configurations A and C consistently decreased as e 
increased. This decrease would be expected because the value of E s decreased 
consistently as e increased (fig. 18). The turns required for recovery for 
configuration B also show an overall decrease, but the trend was not as con- 
sistent as it was for configurations A and C- 


Additional Calculations for Configuration B 


Inasmuch as the turns required for recovery for configuration B did not 
show as consistent a trend with e changes as was noted for configurations A 
and C and inasmuch as configuration B was the only configuration which had an 
oscillatory spin, some additional calculations were made in an attempt to obtain 
a less oscillatory developed spin and to reinvestigate the effects of changes in 
e. In these calculations the estimated values of C m and C? shown in fig- 

r 

ure 20 were used, with all other aerodynamic data inputs unchanged. 


The results of the additional developed-spin calculations for configura- 
tion B are presented in the time histories of figure 21 and in table V. For 
e = 0°, an oscillatory developed spin is still obtained, but the magnitude of 
the oscillations is greatly reduced (compare fig. 13 (a) with fig. 21(a)). 


The results shown in table V indicate that the oscillatory nature of the 
developed spin shown in figure 2l(a) remained about the same for all values of 
e. The average angle of attack remained at approximately 67 ° and the rate of 
rotation decreased as e increased. A summary of these additional results for 
configuration B is presented in figure 22, in which a^verage^ ^e , average; E s , 

and turns for recovery are plotted against e. The magnitude of the spin- 
energy factor E s again decreased consistently as e increased. (The value 
of E s for € = 5° is only 85 percent of the value of E s for e = 0°.) The 
turns for recovery also show an overall decrease as e increases, with a more 
definite trend than that previously obtained when the developed spins were more 
oscillatory. (Compare turns for recovery for configuration B in figs. 1) 
and 22.) Therefore, it appears that the inconsistencies noted in figure 19 in 
the trends of the turns required for recovery for configuration B were caused 
by the more oscillatory nature of the spins from which those recoveries were 
obtained. 
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General Remarks 


The magnitudes of improvement in the spin recoveries shown in tables IV 
and V probably can be achieved with similar configurations having mass distri- 
butions somewhat similar to those used in the present investigation. Any 
attempt to obtain these magnitudes of improvement by readjusting interior mass 
arrangements so as to achieve a larger angle between the XQ-axis (principal 
axis) and the reference X-axis (body axis), particularly angles as large as 5°j 
would, probably be impractical. However, the results presented herein do indi- 
cate that, if a choice exists, the mass distribution which has the largest 
product of inertia will have the best developed- spin and spin-recovery 

characteristics. In addition, the mass distribution with the largest product 
of inertia will have a more oscillatory spin-entry motion which, in some 
instances, may lessen the tendency of a given configuration to enter a spin. 

This trend occurs even though the mass loading parameters which are customarily 
used to predict effects on spins and recoveries remain essentially unchanged. 

In addition, studies made in the past to analyze the spin-entry, developed- 
spin, and spin-recovery characteristics of specific aircraft (refs. 3 and 8) 
have neglected the product-of-inertia term in the equations of motion. As indi- 
cated by the results presented herein, the product of inertia may have an impor- 
tant effect on spins and recoveries and it should therefore be considered in 
future analytical studies pertinent to investigating the spin-entry, developed- 
spin, and spin-recovery characteristics of specific aircraft. 


CONCLUSIONS 


An analytical study of the effects of variations in product of inertia on 
the spin-entry, developed-spin, and spin-recovery characteristics for a delta- 
wing fighter, a swept-wing fighter, and a stub-wing research vehicle with mass 
distributed relatively heavily along the fuselage indicates the following 
conclusions : 

1. For a given mass-loading arrangement such as is customarily described 

in terms of moments of inertia referenced to the body axes, a product of inertia 
due to the angular displacement of the XQ-axis (principal axis) below the ref- 
erence X-axis (body axis) may have important effects on the spin-entry, 
developed- spin, and spin-recovery motions and should be considered in investi- 
gations pertinent to analyzing these motions for specific aircraft. 

2. In general, an increase in product of inertia appears to increase the 
degree of the oscillations occurring during a spin-entry motion. These 
increased oscillations in some instances may lessen the tendency of a given con- 
figuration to enter a spin. For example, for the delta -wing fighter and the 
swept-wing fighter, increasing the angular displacement of the XQ-axis by 5° 
changed the spin-entry motion from a spin entry to a no-spin condition. 

3- Increasing product of inertia generally did not alter the steady or 
oscillatory nature of the developed spins. Increasing product of inertia. 


however, did cause the developed spins to have slightly higher average angles 
of attack, lower rates of rotation, and lower magnitudes of spin-energy factor. 
Because of the lower spin-energy factors, fewer turns were required for spin 
recovery. 

4. As the developed spins become more oscillatory, the effects of 
increasing product of inertia tend to become obscured. 


Langley Research Center, 

Rational Aeronautics and Space Administration, 

Langley Station, Hampton, Va., January 7 7 1965- 



APPENDIX 


EQUATIONS OF MOTION AND ASSOCIATED FORMULAS 


The following equations of motion were used in calculating the spinning 
motions: 


. Iy - Iz , I XZ/* , s pV R 2 Sb 
p = — qr + + pq) + — - 

T X lx 2I X 


+ C J5a 8 » + * ^(c !p P + C !r r) 


T Z - x x 


pr + _ (r 2 . p2) + [C„ + Cmg e & e - £-( 0 ^ ♦ 


pV R Sc 


h_LlX pq + 5£(p - qr) + EXSJI 


Vtd oD -h ( 

C np P + Cn 6a &a + Cn 5r &r + ^(Cn p P + c n r r 


u = -g sin 0 e + vr - wq + 


pV R 2 S 


(°X + c Xs e Se) 


v = g cos 0 e sin 0 e + wp - ur + 


pV R 2 S 


( c YpP + °Y 5a 5 a + C YBr &r) 


• A pV^D S / \ 

w = g cos 00 cos 0e + uq - vr + — ^ — fCz, + 8 e j 


In addition, the following auxiliary formulas were used: 


a = tan"*! — 


p = sin-1 — 


V-p — u 2 + v 2 + w 2 



APPENDIX 


h = u sin 9 e - v cos 0 e sin 0 e - w cos 0 e cos 0 e 


h = hQ + I h dt 


0 e = q. cos 0 e - r sin 0 e 


0 e = p + r tan 0 e cos 0 e + q. tan 0 e sin 0 e 


; 0 e - P 

= — : 

sin 0 e 


0 = sin“-*-^sin 0 e cos 0 e j 


i|r e dt 


Turns in spin = 


2n 
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T ABLE I.- MASS AND DIMENSIONAL CHARACTERISTICS 



C onf i gurat i on 

C onf i gurat i on 

Configuration 


A 

B 

c 

W, lb 

24 8 ll 

23 771 

12 575 

S, sq ft 

695-05 

385.53 

200 

b, ft 

38.12 

35.67 

22.36 

c, ft 

23.755 

11.83 

10.27 

Center of gravity, percent c 

30.0 

55.0 

19.5 

■'•Xq ’ slu S-ft 2 

13 449 

11 709 

4 288 

I Yo , slug-ft 2 

128 000 

82 654 

75 584 

I Zq , slug-ft 2 

138 151 

89 237 

74 867 

1 

H 

-1 022 X 10 _1+ 

-825 X 10 '^ 

-5 556 x 10 -4 

OJ 

'i 



5 e , deg 

-25 

-50 

-50 

S r , deg 

±25 

±6 

± 7 - 

2 

8 a , deg 

±7 

±15 

± 7 - 

'2 
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TABLE III.- MOMENTS AND PRODUCT OF INERTIA 


e , deg 

Tyr, slug-ft 2 

Iy, slug-ft 2 

I z , slug-ft 2 

slug-ft 2 

H C 

1 




mb 



Configuration A 



0 

13 449 

128 000 

138 151 

0 

-1 022 x lO - ^ 

1 

13 487 

128 000 

138 113 

2 176 

-1 022 

2 

13 6oi 

128 000 

137 999 

4 350 

-1 021 

3 

13 791 

128 000 

137 809 

6 518 

-1 019 

4 

Ik 05 6 

128 000 

137 544 

8 677 

-1 017 

5 

14 396 

128 000 

137 204 

10 827 

-1 oi 4 



Configuration B 



0 

11 709 

82 654 

89 237 

0 

-755 x 10 

1 

ll 733 

82 654 

89 213 

l 353 

- 75 ^ 

2 

ll 8 o 4 

82 654 

89 142 

2 704 

- 75 ^ 

3 

ll 921 

82 654 

89 025 

4 052 

-752 

4 

12 086 

82 654 

88 860 

5 395 

-751 

5 

12 298 

82 654 

88 648 

6 731 

-748 



Configuration C 



0 

4 288 

73 384 

74 867 

0 

-3 536 x 10- 4 

1 

4 310 

73 384 

74 846 

1 232 

-3 53 ^ 

2 

4 374 

73 384 

74 781 

2 462 

-3 531 

3 

4 481 

73 384 

74 674 

3 689 

-3 526 

4 

4 631 

73 384 

74 524 

4 911 

-3 518 

5 

4 824 

73 384 

74 331 

6 128 

-3 508 
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TABLE IV.- RESULTS OF DEVELOPED -SPIN CALCULATIONS 


e, deg 

a, deg 

deg 

rad/sec 

Vr, ft/sec 

E s 

Turns for 
recovery 


Configuration A 


73.8 

-2.0 

1.26 

318 

0.1308 

3.44 

74.0 

-1.9 

1.24 

318 

.1257 

3-33 

74.0 

-1.9 

1.21 

318 

.1185 

3.21 

74.1 

-1.8 

1.19 

318 

.1136 

3.10 

74.1 

-1.8 

1.17 

317 

•1093 

3-02 

74.1 

-1.8 

1.15 

317 

.1044 

2.93 


Configuration B* 


0 

65.6 to 85.7 

-25.3 to 17.6 

1.69 to 2.63 

341 to 343 

0.4213 

5.66 

1 

63.7 to 87.3 

- 25*5 to 17.3 

1.59 to 2.68 

344 to 346 

.4007 

4.10 

2 

6l.O to 88.5 

-17.I to 13.8 

1.58 to 2.62 

345 to 348 

.3785 

4.37 

3 

60.5 to 90.7 

-l 6.4 to 15*0 

1.49 to 2.6l 

345 to 348 

.3601 

3.36 

4 

61.9 to 92.8 

- 15*1 to 11.9 

1.45 to 2.57 

345 to 347 

• 3513 

3.44 

5 

59*8 to 92.5 

- 14.7 to 14.5 

1.39 to 2.62 

345 to 348 

*3396 

3*54 



Configuration C* 




0 

83.6 to 87.2 

- 6.4 to 0.2 

2.07 

277 

1.5967 

6.98 

1 

84.5 to 87.3 

-5-6 to .1 

1.98 

276 

1.4694 

6.22 

2 

85.1 to 87.4 

-4.9 to -.3 

1.91 

276 

1.3646 

5.88 

3 

85.7 to 87.8 

-4.3 to -.3 

1.84 

276 

1.2642 

5.66 

4 

85.1 to 88.2 

- 4.0 to -.4 

1.74 

276 

1.1256 

5.28 

5 

85.5 to 88.6 

-3.5 to -.4 

1.65 

275 

1.0169 

5-01 


*Range of oscillations given. 
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TABLE V.- RESULTS OF ADDITIONAL DEVELOPED-SPIN CALCULATIONS 

FOR CONFIGURATION B* 


e , deg 

a , deg 

( 3 , deg 

rad/sec 

Vp, ft/sec 

E s 

Turns for 
recovery 

0 

. 

61.5 to 71.6 

- 5.9 to 0.3 

1.11 to 1 . 4 l 

359 to 362 

0.1175 

3.67 

1 

62.4 to 71.8 

- 4.7 to 0 

1.11 to I.37 

359 to 361 

.1135 

3.64 

2 

62.6 to 72.5 

-5.2 to 0.3 

1.08 to I.33 

359 to 361 

.1064 

2.92 

3 

62.0 to 73.3 

- 4.9 to 0.4 

1.05 to I.35 

! 358 to 361 

. 1045 

2.90 

4 

62.2 to 74.0 

-5-5 to 0.8 

1.02 to I.35 

357 to 360 

.1017 

2.26 

5 

62.6 to 74.8 

- 5.5 to 1.1 

1.00 to I.35 

357 to 360 

• 0995 

2.29 


-X- 

Range of oscillations given. 
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Projection of 
relative wind 



(b) 0 e and. 0 e = 0. 



(c) 0 e and i|r e = 0, and in this case 0 = 0 e . 

Figure 1.- Three-view sketch showing body system of axes, related angles, and positive 
directions of corresponding forces and moments. 
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Figure 3-- Variation of static longitudinal stability characteristics. Measurements are 
referred to 30 - 9 33-^ and 19 -percent mean -aerodynamic -chord positions for configurations 
A, B, and C, respectively. 





e 4. - Variation of static lateral and directional stability derivatives 
tack. Measurements are referred to the 30- 33- > and. 19-percent mean-a' 
sitions for configurations A, B, and C, respectively. 










<$> e , deg 20 



Figure 12.- Calculated developed-spin and spin-recovery motions for configuration A 





(a) € = 0°; Ixg = 0 slug-ft 2 . 

Figure 13.- Calculated developed- spin and spin-recovery motions for configuration B 





(b) e = 5°; Ixz = 6 731 slug-ft 2 . 


Figure 13 . - Concluded. 





(a) e = 0°> 1 x 2 , = 0 slug-ft 2 . 

Figure 14.- Calculated developed- spin and spin-recovery motions for configuration C 





average 

sec 


Configuration 



e } deg 


Figure l6.- Variation of developed-spin average 


rate of rotation with e for all configurations 







4 =* 

4 =- 



0 r 



CL ; deg 


Figure 20.- Estimated variation of C m 


and C 7 coefficients with 
figuration B. 5e = -20°. 


a, used for additional calculations for con- 


[ 






70 


60 


0 I 


1.2 


^e, average' g 
rad/sec 


.4 


0 


E 


s 


.4 


0 



Turns 

for 

recovery 


2 




1 h 


0 


1 


1 

2 


3 


4 


€, deg 


Figure 22.- Variation of developed- spin average angle of attack, average rate of rotation, 
spin-energy factor, and turns required for recovery with e for additional calculations 
for configuration B. 
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"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenome?ia in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof." 

— National Aeronautics and Space Act of 1958 
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